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Abstract 

Normal state orbital diamagnetism of charged bosons quantitatively accounts for recent high- 
resolution magnetometery results near and above the resistive critical temperature T c of supercon- 
ducting cuprates. Our parameter-free descriptions of normal state diamagnetism, T c , upper critical 
fields and specific heat anomalies unambiguously support the 3D Bose-Einstein condensation of 
preformed real-space pairs with zero off-diagonal order parameter above T c at variance with phase 
fluctuation scenarios of cuprates. 
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A possibility of real-space pairing, as opposed to the Cooper pairing, has been the subject 
of many discussions, particularly heated over the last 20 years after the discovery of high 
temperature superconductivity in cuprates [l|. The first proposal for high temperature su- 
perconductivity, made by Ogg Jr in 1946 Q], already involved real-space pairing of individual 
electrons into bosonic molecules with zero total spin. This idea was further developed as 
a natural explanation of conventional superconductivity by Schafroth and Butler and Blatt 
P). However, with one or two exceptions, the Ogg-Schafroth picture was condemned and 
practically forgotten because it neither accounted quantitatively for the critical behavior 
of conventional (i.e. low T c ) superconductors, nor did it explain the microscopic nature of 
attractive forces which could overcome the Coulomb repulsion between two electrons con- 
stituting a pair. The failure of the 'bosonic' picture of individual electron pairs became 
fully transparent when Bardeen, Cooper and Schrieffer ^| proposed that two electrons in a 
superconductor were indeed correlated, but on a very large distance of about 10 3 times of 
the average inter-electron spacing. 

Highly successful for low-Tc metals and alloys the BCS theory has led many researchers to 
believe that novel high-temperature superconductors should also be "BCS-like". However, 
the Ogg-Schafroth and the BCS descriptions are actually two opposite extremes of the 
same electron-phonon interaction. Indeed by extending the BCS theory towards the strong 
interaction between electrons and ion vibrations, a charged Bose gas (CBG) of tightly bound 
electron pairs surrounded by lattice deformations (i.e. of small bipolarons) was predicted 
by us with a further prediction that high T c should exist in the crossover region of 
;he electron- lattice interaction strength from the BCS-like to bipolaronic superconductivity 
f|. Experimental evidence for an exceptionally strong electron-phonon interaction in novel 
superconductors tIIs!, 3.1lo| is so overwhelming that bipolaronic CBG [ill] could be a feasible 



alternative to BCS-like scenarios of cuprates. Nevertheless, some authors [12] have dismissed 
any real-space pairing, advocating a collective pairing (i.e Cooper pairs in the momentum 
space) at some high temperature T* which are phase coherent below a lower temperature 
T c < T*. Ref.j3| has argued that the superconducting transition in cuprates is an almost 
two-dimensional Kosterlitz-Thouless (KT) transition, where a vortex liquid exists above T c 
different from the BCS theory and its strong-coupling bipolaronic extension with a perfectly 
"normal" state without any off-diagonal order. 

So far there has been no decisive conclusion on the origin of anomalous normal state 
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of cuprates. Some normal state properties have been satisfactorily interpreted within the 
7ermi-liquid approach, while many others have been understood with preformed real-space 
111 ] or Cooper [3| pairs, in particular on the underdoped side of the phase diagram. More- 
over preformed real-space pairs could coexist with the Fermi-liquid, which effectively hides 
them in the normal state kinetics. Any direct evidence in favor of either scenario is highly 
desirable. If real-space pairs indeed exist in superconducting cuprates, then their supercon- 



ducting state should 
Its critical behavior 



je a three-dimensional (3D) Bose-Einstein condensate (BEC) of CBG. 



BCS 



ll| is rather different from any "universal" criticality like mean- field 



3D "XY" or KT 



12] transitions. 

Here I show that high-resolution magnetometery in the critical and normal regions pro- 
vides unambiguous evidence for real-space pai ring in cuprates. 



A number of experiments (see, for example, 
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121 Q 

and references therein) , 



including torque magnetometery, showed enhanced diamagnetism above T c . Originally it was 
explained as the conventional fluctuation diamagnetism in quasi-2D BCS superconductors 
(see, for example Ref. Q|). The data taken at relatively low magnetic fields (typically below 
5 Tesla) revealed a crossing point in the magnetization M(T, B) of most anisotropic cuprates 
(e.g. Bi-2212), or in M(T, B)/B^ 2 of less anisotropic YBa 2 Cu 3 7 _ 6 Q- The dependence of 
magnetization (or M/B 1 / 2 ) on the magnetic field was shown to vanish at some characteristic 
temperature below T c . Importantly more recent data taken in high magnetic fields (up to 
30 Tesla) show that the crossing point, anticipated for low-dimensional superconductors 
and associated with conventional superconducting fluctuations, does not explicitly exist in 
magnetic fields above 5 Tesla jl£ , Q . 

Ref.j3| has linked the enhanced normal state diamagnetism with mobile vortexes well 
above T c where conventional fluctuations should be negligible. Surprisingly the same torque 



magnetometery 
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ll^ uncovered that the diamagnetic signal above T c increases in 
magnitude with applied magnetic field, B. Such magnetic field dependence of magnetisation 
M(T, B) is entirely inconsistent with what one expects from vortex liquid. While —M(B) 
decreases logarithmically at temperatures well below T c , the experimental curves clearly show 
that —M(B) increases with the field at and above T c , just opposite to what one could expect 
in a conventional vortex liquid. These significant departures from the London liquid behavior 
indicates that vortex liquid does not appear above the resistive phase transition (see also 
Ref . [l^ ) . Also accepting the vortex scenario and fitting the magnetization data in Bi-2212 
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FIG. 1: Bipolaron picture of high temperature superconductors. A corresponds to a singlet oxygen 
intersite bipolaron, B is a triplet intersite bipolaron. 



with the conventional logarithmic field dependence |18f| , one obtains surprisingly high upper 
critical fields H c2 > 120 Tesla and a very large Ginzburg-Landau parameter, k — Ah/£ > 450 
even at temperatures close to T c . The in-plane low-temperature magnetic field penetration 
depth is \h ~ 220 nm in optimally doped Bi-2212 (see, for example Hence the zero 

temperature coherence length £ turns out to be about the lattice constant, £ < 0.5nm. 



Such a small coherence length is per 



ectly compatible with direct STM measurements of the 



individual vortex cores in Bi-2212 3] and with the size of the vortex core in CBG [ll | 
However it rules out the "preformed Cooper pairs" |l2j, since the pairs are virtually not 
overlapped at any size of the Fermi surface. 

Here I calculate the magnetization, M(T,B), of anisotropic CBG on a lattice, and com- 
pare the result with diamagnetism of cuprates recently measured in Ref . [l8( . A low-energy 
structure of cuprates in the bipolaron model is shown in Fig.l, where oxygen holes are 
bound into real-space intersite singlets (A) and triplets (B) separated by an exchange en- 
ergy J jscj, which is estimated as a few tens or hundreds Kelvin depending on doping 



in agreement with experimental charge and spin pseudogaps in cuprates [2lj . Bipolarons 
are almost ideal charged bosons, because their Coulomb repulsion is strongly suppressed 
by a large lattice dielectric constant. When the magnetic field is applied perpendicular to 
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copper-oxygen plains the quasi-2D bipolaron energy spectrum is quantized as 

E a =w(n+ 1/2) + 2t e [l-coa(K z d)], (1) 

where a comprises n = 0, 1, 2, ... and in-plane K x and out-of-plane K z center-of-mass quasi- 
momenta, uo = 2fteB / \jm x m y) t c and d are the hopping integral and the lattice period 
perpendicular to the planes. The spectrum consists of two degenerate brunches, the so- 
called "x" and "y" bipolarons [llj . with anisotropic in-plane bipolaron masses m x = m and 
m y m 4m. Expanding the Bose-Einstein distribution function in powers of exp [(fj, — E) /k B T] 
with the negative chemical potential fi one can after summation over n readily obtain the 
boson density 

^ _2eB^ T ln+ _ /u ^exp[(/i - 2t c )r/k B T] 



^EVW)=^, (2) 
7r "" ~l 1 ~ exp(— tur/ k B T) 



and the magnetization, M(T,B) = — k B Td/dB J2 a m [l — exp(yU — E a )/kBT], 



2ek B T ^ T ( 2t c r 



r=l 



M(T,B) = -n m+ ^^I^)x (3) 



exp[(/i — 2t c )r/k B T] (\ lo exp(— ur/k B T) 



1 — exp(— ur/k B T) \r k B T[l — exp(—ur /k B T)] 
where fib = fte/ ^Jm x m y , ft — fi — ou/2 and Iq(x) is the modified Bessel function. At low tem- 
peratures T — > Schafroth's result jsj is recovered, M(0, B) = —Ubfib- The magnetization 
of charged bosons is field-independent at low temperatures. At high temperatures, T>r c 
the chemical potential has a large magnitude, and we can keep only the terms with r = 1 in 
Eqs.(2,3) to obtain M(T,B) = —nbHbU)/(Gk B T) at k B T ^> k B T c ^> lu, which is the familiar 
Landau orbital diamagnetism of nondegenerate carriers. 

The critical region r = T/T c — 1 < 1 requires numerical calculations, which have been 
done H for an anisotropic 3D CBG with t c > k B T c /2 and I Q (x) « e x /\[2^x in Eqs.(2,3). 
Notwithstanding, one can nicely map the exact results, Fig.2, with a simple analytical 

expression by replacing summation over all but the first Landau level for integration, 

-l 

(4) 



M(T, B) 0.22u 



0.37cj 
k B T c 



r+ W— — + r 2 



n b fi b k B T c 

Remarkably Eq.(4) predicts almost field-independent diamagnetism well below T c , |r| ^> 
u/k B T c , a linear field dependence M(T C , B) ~ B well above T c , r 3> u/k B T c , and an unusual 
square root behavior at T = T c , M(T C ,B) ~ B l l 2 . Here T c is the familiar Bose-Einstein 
condensation temperature k B T c = 3.31ft? (nb/2) 2 ^ 3 / ' (m x m y m c Y^ 3 ', with m c = ft? /\t c \d 2 . 
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FIG. 2: Exact numerical magnetization (symbols) compared with the analytical approximation, 
Eq.(4) (lines) . 
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FIG. 3: Diamagnetism of optimally doped Bi-2212 (symbols) 18] compared with magnetization of 
CBG, Eq.(6), near and above T c (lines). 



Comparing with experimental data one has to take into account a temperature and field 
depletion of singlets due to their thermal excitations into spin-split triplets and single polaron 
states, Fig. IB. If the spin gap is small compared with the charge pseudogap, J < A/2, 
triplets mainly contribute to temperature and field dependencies of the singlet bipolaron 
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density near T c , 

n b (T,B) = n 6 (T c , 0)[1 -or- (B/B*)% (5) 

where a = Z(2n c t)- 1 [J(e J / kBT - - l)" 1 - k B T c \n(l - e - J / fcsTc )], fi B B* 
(2k B T c n c t) 1/2 sinh(J/ '2k B T C ), \i B m 0.93 x 10~ 23 Am 2 is the Bohr magneton, n c is 
the atomic density of singlets at T = T c in zero field (n c < 0.1 in optimally doped 
cuprates), and 2t is the triplet bandwidth, which is taken much larger than k B T c . A triplet 
contribution to diamagnetism remains negligible compared with the singlet diamagnetism 
if r J/k B T c . Then Eq.(4) mapping numerical magnetization in the critical region is 
modified as 



, (B/B*) 2 I 0.37B/B f (B/B*) 2 \ 
1 + 2a/3 V (1 + 2a/3) 2 V 1 + 2a/3 J 

where B = k B T c /2jj, b . Using the magnetic field in-plane penetration depth, A^ 2 ~ 21 



; magnetic neia m-pi 
Q and of CBG H, 



(fim) 2 of optimally doped Bi-2212 j23j and of CBG \^ = 2nbe 2 (rn x + m y )/([i 'm x my), 
we estimate the bipolaron mass asm^ 7.5m e in agreement with the analytical and numerical 
QMC results H, and n 6 (T c ,0)/i b = h^o(m x m y ) 1 / 2 /2e\ 2 H (m x + m y ) w 2100A/m. The BEC 
temperature corresponds to the temperature were the in-plane resistivity starts to drop 
with temperature lowering, which is about T c = 90K in optimally doped Bi-2212, so that 



18], 



B = 524 Tesla. This choice of T c = 90K is also justified by low-field magnetization 
which has an exponent close to 1/2, M(90K, B) ~ B 1 / 2 just for this temperature. The 
remaining two parameters in Eq.(6) are found using the experimental field dependence of 
M(T, B) at any fixed temperature near T c . Fitting M(T, B) at T = 89 K, Fig.3, yields 
a = 0.62 and B* = 56 Tesla, which according to Eq.(5) corresponds to the singlet-triplet 
exchange energy J w 20K. Quite remarkably all other experimental curves in the critical 
region are well described by Eq.(6) without any fitting parameters, Fig.3. 

I conclude that the normal state diamagnetism observed in many cuprates [3, Q, [3, 
llil IttI [3] provides unambiguous evidence for charged real-space bosons. The exper- 
imental data, Fig.3, clearly contradict BCS (with or without conventional fluctuations) 
and KT scenarios of the phase transition in cuprates. If we define a critical exponent as 
5 = In B / In \M(T, B) \ for B — > 0, the T dependence of S(T) in CBG is dramatically different 
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rom KT and other "universal" critical exponents, but it is very close to the experimental 
3 S(T), Fig.4. 

Another strong argument in favor of 3D BEC in cuprates has been drawn using parameter- 
free fitting of experimental T c with BEC T c in more than 30 underdoped, optimally and 
overdoped samples |24|. Whereas the KT critical temperature expressed through the in- 
plane penetration depth Q ^bTkt ~ 0.9dh 2 /(16Tce 2 \j { ) appears several times higher 
than the experimental values in many cases. There are also quite a few samples with about 
the same Xh and the same d, but with very different values of T c , in disagreement with the 
KT transition. The large Nernst signal, allegedly supporting vortex liquid in the normal 



state of cuprates 



261 ] . has been explained as perfectly normal state phenomenon owing to a 



partial localization of charge carriers in a random potential inevitable in cuprates |27j |. CBG 
upper critical field and the specific heat in the magnetic field have been found in striking 
consensus with experimental data j3] following our prediction [29]. More recently the d- 
wave symmetry and real-space modulations of the order parameter have been also explained 
with CBG in underdoped [ll] and overdoped 30] cuprates. 

I thank Victor Kabanov for his help with numerical calculations and valuable discussions. 
The work was supported by EPSRC (UK) (grant EP/C518365/1). 
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